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The groimdwave pulse has been considered in detail by Johler [1962, 1963b]. A theory 
was presented as a direct numerical and analytic evaluation ot the Fourier transform in- 
tegral. Such a theory can be extended to treat pulsed ionospheric waves. 

The reflection process for various model anistoropic ionospheres has been studied in 
detail in previous papers by Johler [1962]. The analysis of the reflection coefficient thus 
demonstrated is employed to construct transfer functions for the ionospheric wave. These 
transfer functions are then introduced into the Fourier transform of the pulse, and the 
behavior of LF pulses propagated by a reflection process at the ionosphere is determined for 
three models for the lower ionosphere. The composite ground and ionospheric pulse is 
formed, and a method of tagging a point-in- time on the pulse is detailed. 



1. Introduction 

The low frequency radio signal [Johler, 1962] propagated from a transmitter at the surface 
of the terrestrial sphere can overcome the curvature of the earth and reach a receiver some 
distance away on the surface of the earth by diffraction as a groimdwave pulse or by a single 
reflection or multiple reflections at the lower ionosphere. Numerous investigators have treated 
the groimdwave pulse or signal theoretically [Novikov, 1960; Johler and Walters, 1959; Keilson 
and Row, 1959; Pekeris and Alterman, 1957; Wait, 1955 and 1959]. Thus, for example, 
Wait and Froese [1955]; Wait and Conda [1959] have developed closed form solutions for 
certain pulses for relatively simple situations. On the other hand, very little has been accom- 
plished on the pulse signal reflected at the ionosphere at LF. A procedure was outlined by 
Johler [1962] for treating theoretically both the groundwave and the ionospheric wave pulses 
by a direct numerical integration of the Fourier transform-integral. Thus, regarding the 

[5 -> -» 5 -» -» 

continuous time-harmonic-wave solution of Maxwell's equations, -r- E=ic*)E: — H=io)H. as 

1 ot dt 

the transform or transfer function of the propagation medium, it is no longer necessary for the 
case of linear amplitude transforms to resolve Maxwell's equations for each pulse. It is only 
necessary to form the product of the pulse transform and the complex transfer function of 
the medium and integrate to the time domain. 

The reflection process for various model anisotropic ionospheres has been studied in con- 
siderable detail by Johler [1962, 1963a; Johler and Harper, 1962a and 1962b]. Results obtained 
with the aid of this analysis are employed in this paper to construct three transfer functions 
for three models of the propagation medium. These transforms are then introduced into the 
integrand of the Fourier integral of specified pulses, and the behavior of such pulses propagated 
by reflection at the lower ionosphere is demonstrated. The composite ground and ionospheric 
pulse is formed, and the methods for tagging a point-in-time [Johler, 1963b] on the pulse are 
detailed. 



i This work was sponsored by the Advanced Research Projects Agency under ARPA Order No. 183-02, Amendment 4, or NBS Project 85411— 
initiated under this contract. This work forms the theoretical basis for the experimental program pursuant under the same contract on NBS 
Project 85461. 
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2. Theory of Transient Propagation 

A radio signal, 2?(£, d) } v/m or i?(£, <i) amp-turns/m, at a time, £, and a distance, r/, along 

the surface of the terrestrial sphere relative to a transmitter is related to the continuous time 

-> -> 

harmonic solution of Maxwell's equations, E(a), d) or H(co, d) at a frequency /=co/27r c/s by the 

Fourier transform-integral for linear amplitude restrictions, 



-j /» 00 /» 00 

E{t',d)=j- exp (iut')E(w,d)f r (a) exp (-iut)F s (t)dtda 



(1) 



where / r (co) is the transfer function of the receiver or measuring device and 2^(0 is the transient 
dipole source current moment, and the local time, t', is defined, 



t' = t — r)id/c, 



(2) 



where 171 is the index of refraction of air (77^ 1.000338 at the surface of the earth), and c is the 
speed of light (c^2.997925(10 8 ) m/sec). The propagation medium transfer function, E(w, d), 
can be split geometrically-optically [Johler, 1962] into the sum of a multiplicity of ordered 
rays, j=0, 1, 2, 3 . . ., 

tf(«,dHS^(«,4), (3) 

where j=0 corresponds to the ground wave, j=l corresponds to the first ionospheric reflection, 
etc. As a consequence of these rays, the pulse fields exhibit a corresponding series, 



E(t',d)=y: l E j (t' > d), 

}-0 



(4) 



separated in time by ever later times, V ', as the order, j, is increased. The groundwave, j=0, 
and the first ionospheric reflection will be considered in this paper. Obviously the latter 
(j=l) will arrive at the receiver later, since the ray length between the transmitter, the iono- 
sphere, and the receiver is longer than the distance, d, along the surface of the earth, to the 
receiver appropriate to the groundwave, j=0. These facts have been verified experimentally 
[Johler, 1962], and the ray theory for propagation between two points, s (source), and (ob- 
server) on the surface of the earth can be summarized as shown diagrammatically, figure 1. 
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Figure 1. Propagation media and boundaries, illustrating various geometric-optical rays in the vicinity of the 

terrestrial sphere. 
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The electrical properties of the three principal media of propagation are described by the wave 
numbers, k h k 2 , k h 



CO 



ki=- 7/i (5) 

c 



C \ co 

£3=" utf (7) 

where e 2 is the dielectric constant relative to a vacuum, e 2 eo=e the permittivity, e =-^— ■> and 

Mo=47r(10 -7 ) h/m, the permeability of space; a, mhos/meters is the conductivity of the 
ground. The index of refraction for the ionosphere, r\l\ r e , is quite complicated, and there 
are four such indices for a homogeneous layer, ordinary (0) and ext inordinary (e), upgoing (i) 
and ordinary (0) and extraordinary (e), downgoing (r) indexes. Furthermore, the electrical 
properties of k 3 vary with altitude, h,k 2 = k 3 (h) which fact is taken into account in the reflection 
coefficient, T, of the ionosphere. The lowest electron density considered is taken as the lower 
boundary at the altitude, h, for purposes of geometric-optical rays, figure 1. 

The groundwave (j=Q) signal, E (t' , d) is considered to be the pulse which travels from 
s to 0, figure 1, along the surface of the earth taking into account the electrical properties of 
the earth, k 2 , and the medium, k x , assumed to be of infinite extent. This is the earliest signal 
to arrive at 0, and the earliest time such a signal can arrive at is £'=0 or t = rjid/c. Somewhat 
later (about 20 to 40 /*sec) the signal propagated over the ray, i=l, and reflected from the 
lower ionosphere at a height, h, and angle of incidence, </> u , will arrive at a local time, £j=£i^>0, 
or will take on appreciable amplitude at a time, £j>0. Obviously, higher order rays 0=2, 
3, 4 . . .) are possible and indeed have been observed [Johler, 1962]. However, for purposes 
of illustrating technique the distances, d=S00 to 2000 statute miles or [1200 to 3000 km will 
be considered in which for the assumed models, j=Q and j=l, are of prime consideration. 
Indeed, a model will be considered in which the j=l jay, figure 1, is diffracted beyond the 
geometrical-optical horizon, or beyond the point at which the ray is grazing incidence on the 

earth, r j = r 1 =-' The rays reflecting to and fro between the earth and the ionosphere are 

L 

influenced by the local reflection coefficient at the ground surface, r=a, and at the ionosphere 
surface, r=a-\-h. Distance, d = ad, increases in the positive ^-direction or the signals are as- 
sumed to be progressing toward the receiver or observer, 0, in the positive ^-direction, figure 1. 
The transform for the propagation medium, considering the groundwave, Eq(q) } d) } and the 
first ionospheric wave, Ei(w,d) can be written [Johler, 1962] for vertical dipoles at s 
and [Johler, 1961]. 

E(a,,d)==E (o>,d)+E 1 (a,,d) (8) 

,,d)=K«C W'Ua) 1 ' 8 ^ !- -* ^ " 2a 4JJ - 

L a J V. [AT s — i/b e \ 

F\T ee I 



E(c 

l-iadD^C exp [i^-t'WiaJ^Tce V (9) 



where 

C=(10- 7 )Iol/d 

and 

(7^ sin n, 



r = t- Vl DJc, (10) 

Di=2[(a-\-h) cos fa—a cos n], (11) 

fa=fa,i, (12) 
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r=n. (13) 

a~0.85, the earth radius factor, a is the radius of the earth, t s are the special roots of Riccati's 
differential equation [Johler, "Walters, and Lilley, 1959], 

^-25 2 r+l = (14) 



°e~ rU =TT7i' CIS; 



where for vertical polarization 8 = 5 

(^[§-1] 

where h and k 2 are defined in (5) and (6) , 

. d 
a sin - 

sin<p«== — spy-, (16) 

2a(a+h) (l-cos -)+h 2 

=(l+^) [2 sin ^sin dj n JTa (l-cos £)+/f]/[(a+^) cos |- a ] V'^!, (17) 

^~(| z^' 2 HZ{z x ) exp { -tIfcr/12-*] } (18) 

S!= fe a cos 3 r0/(3 sin 2 r,) (19) 

J «, exp [-«2W3] yy l {p)—gWi{p) 



«1 



where 6' = (d—d H )/a, d H is the distance to geometric-optical horizon, figure 1, ri=~; 



2= -^TJ ^V 1 -^ (21) 

[Wait and Conda, 1958; Johler, 1962]. Wi(z) and TFi(g) are airy integrals related to the 
Hankel functions, Hl%(z), m%(z), 

W 1 (p)=exp [-27ri/3]V^/3(-p) I/2 ^|[!(-p) 3/2 ], (22) 

Wi(p)=exp [-4™/3]y£/3ptf^[f(-p) 3/2 ], (23) 



where 



7T . .57T 

— g ^ arg p<y , 



and arg (— p)=arg p— t, arg p m/n = — arg p where m and n are integers. F x can also be calcu- 

lb 

lated in the geometric-optical diffraction zone (ti=-> fig. 1) with the aid of [Wait and Conda, 
1958], an alternate formula, 

F^-2i^ exp [-***]§ (X^flffi ( 24 ) 

This formula clearly shows the groundwave nature of the diffracted wave, since it is the classical 
residue series by analogy to first term of (8) and (9), which is the groundwave. Hi 2 f l(z) was 
computed with the aid of a computer technique [Berry, 1962] using the relations, 

#$(*)={ -exp (^) J lfl (2)+J- W («) \ji sin (-|) (25) 

#$(z)={-exp (?£) J w (z)+J- W (z)\ji sin (-^) (26) 
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where 

(—l) m (z/2) v+2m 



Jv{z)=H 



(z/2y r (z/2y 



2" (o+l) 



r _ (s/2)"» (g/2)' (s/2)« "] / 1 2\ 

L l!(t7+l)" r '2!(«7+l)(t7+2) 3l(»+l)W)(»+3) t,,, JV 3' 3/ l ; 



e/t,(2) is made single valued by specifying arg {zj2) v — v arg (z/2) for A— 2r<arg (s/2) <A. Ordi- 
narily, A=ir but may have different values, depending upon applications. 

As \z\ gets large, (27) converges slowly and loses precision. Also, the precision lost in 
(25) and (26) is serious, since in half of the 2-plane the terms on the right are nearly negatives 
of each other when Im (z) is large. It is therefore necessary to use the asymptotic scries, 



where, 



g ( 2 ) (s)= exp [-fr+jg+frrfll Ti(JMg)> ( _ 2?r<ftrg 2<7r) 

■y/TZ/2 

r.(*)=i+g- nl(4g)ll - (28) 



T c (2) is asymptotic except for v=n-\-% and integral n, and hence should be terminated at the 
smallest term, the next term being greater than the error of such a termination. 

The reflection coefficient, T ec , of the ionosphere comprises a quite complicated computa- 
tion which has been detailed by Johler [1962] and Johler and Harper [1962a and 1962b]. The 
ionosphere is considered locally, figure 1, at the reflection point of the ray; that is, a point below 
which the electron density <^5-10 electrons/cm 3 . This point in the ionosphere can be tested 
for contribution to the reflection coefficient by varying the lower boundary. At LF and VLF 
the reflection coefficient will be found to be determined by greater values of electron density. 
A plane wave, E t , is assumed to be transmitted into the plasma, 

E t = \E\ exp [i (««-» nDjj, (29) 

where 

r]D=x sin (j>i sin <£ a +2/ sin <f> t cos <£ a +2f, (30) 

in which f is in general a complex number related to the index of refraction, 77, 

^ 2 =f 2 +sin 2 (/>,, (31) 

and a right-handed x, y, z Cartesian coordinate system is assumed locally at the lowermost 

ionosphere boundary, so that the y-z plane is oriented toward magnetic north and, hence, 

-> 
contains the earth's magnetic field vector, &{ . The angle <£ a , is measured clockwise toward 

the x-z plane from magnetic north as the orientation of the plane of incidence of the incident 

wave, and the dip or inclination angle, /, is measured downward from the horizontal x-y plane. 

The angle of incidence, <fo, is measured from the vertical. 

The quantity, f, is determined by a simultaneous solution of the continuous time-harmonic 

form of Maxwell's equation, 

~ E=ia>E, ^ H=iuH, VXE+n Q ia*H=0, VXH-J-e ia>E=0, (32) 

and the equation of motion of an electron, 



-> -> 



iwmV+mgV+i*>e(yx&()+eE=0, (33) 
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and J= —NeV for N electrons/cubic meter of charge, e, and mass, m. V is the electron velocity 
vector. 

Equations (32) and (33) imply that Maxwell's equations can be written in terms of the 
dielectric tensor, ||e||, 

VXE+» ia>H=Q 



VXH-\\e\\-E=0 
where in Cartesian x, y, z coordinates, Johler and Harper [1962b] 



(34) 



c|| = 



e ^a; 



where 



1+: 



ih L 



— ih 7 



-ih L 



ihi 



s 2 -h 2 

~s(s 2 -h 2 ) 
h L h T 



s(s 2 -h 2 ) 



s 2 -h 2 

s 2 -hl 

s{s 2 -K 2 ) 

s 2 -hl 
s(s 2 -h 2 )_j 



(35) 



<*N 



!['-=} 



i, 



= — 2~> 



h L =—h sin I, h T =h cos 7, o: H =fi e^{ /m, !£( = ]${ \, and oo^=Ne 2 /e 7n, and j , H =o) H /2ir, and 
j n = u n j2tt are the gyro and plasma frequencies, respectively. SK =H m ^0.5 gauss or 40 amp- 
turn s/m. 

The complex, frequency-dependent parameter, g, in the equation of motion (33) is found 
by an integration over a Maxwellian energ} T distribution expressed in terms of energy, eu=^rnV 2 , 
or fo— [ejirkTY' 2 exp [—eu/kT], in which there exists momentum transfer collisions with gas 
molecules proportional to the energy, u, v=v(u) . Here kT is Boltzmann's constant-temperature 
product. Such an integration can be readily accomplished for a continuous electromagnetic 
wave of angular frequency, o)=2irf. Thus, 

K~ 4 f r rsfs] £*{(ir) ' exp (-«^*r)*}] '-* (se) 

where 12 is an angular frequency yet to be identified [Johler, 1962; Johler and Harper, 1962a 
and 1962b]. 

The dielectric tensor can therefore be generalized to take account of collisions, v=v(u), 
proportional to energy [Johler and Harper, 1962b, 



gWiHifii g(Q 2 )+ifi s 



sin I 



j sin I Zl. sin I 

gfOiHiOi g(Q 2 )+iQ 2 



2 [ t- cosI 4- cos 1 



g (Q 1 )+in 1 g(n 2 ) + in £ 



■ COS I — COS I 



iOp + < 



i sin 2 I i sin 2 I 



g(fi )+io g(n^+io x ) g (n 2 )+ic 



g(fii)+ifli g(fi 2 )+iC : 



-sin I cos I 



■ sin I cos I 1 sin I cos I 



g(n ) + ifi g(^) + in, g (n 2 ) + i^ 



- sin I cos I I sin I cos I | sin I cos Ij 
g(fio) + i^o + g(0i) + ifii + g(n 2 ) + i P, 



3 I 1 cos 2 I 



g (n )+io gtfij+ip.! g(n 2 )+in 2 



where three characteristic frequencies, 12, have been identified, 12 = co, Q 1 = o-\-o) H , and 2 =co— «#. 

-> -> 

The vectors, F and iJ, can be eliminated, whereupon it can be concluded a quartic equation 
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in f is satisfied if the. field, E, exists, 

a»f*+arf*+arf , +a 1 r + (J. =0. 



(38) 



Equation (38) is analogous to the Booker [1939] quartic in classical magneto-ionic theory 
and can be readily solved Ity various numerical techniques [Johler and Harper, 1962a and 
1962b]. Since the quartic has four roots, the four indexes of refraction of the ionosphere, 
Tjolo can be identified with the aid of (31). 
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Figure 2. Model ionosphere profiles, N(h), illus- 
trating the variation of the electron density (elec- 
trons/cm z ) of the lower ionosphere with altitude: (a) 
Houston Composite model appropriate to lower geo- 
graphic latitudes; (b) Blackout Model appropriate 
to higher latitudes (with Gaussian distribution in 
the region J^O-60 km) . 
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-ORDINARY (o) MODE 

- EXTRAORDINARY ( e ) MODE 
UPGOING (i) AND DOWNGOING (r) 
WAVES ARE COUPLED AT EACH BOUNDARY 

Figure 3. Structure of the stratification employed in 
the analysis procedure, illustrating the concept of 
continuous stratification. 

The ordinary and extraordinary upgoing and downgoing propaga- 
tion components are coupled at each boundary. 



The detailed structure of two model ionospheres is illustrated, figure 2, as the electron 
density-altitude, N(h), profile. The collision frequency-altitude, v(h) profile has been previously 
described as the Nicolet/3 profile, Johler [1962]. These structures were represented by a 
flexible plasma stack of slabs illustrated, figure 3 of arbitrary thickness (except for the top- 
most slab of thickness, z p =co). Both the number of such slabs, p, and the thickness were 
quite flexible, since the idea of this analysis is that the measured electron density-altitude and 
collision frequency-altitude (N(z) and v(z), respectively) profiles can be approximated to any 
desired accuracy by decreasing z n and increasing p simultaneously until a stable reflection 
process is obtained. 

A constant electron density, collision frequency and static magnetic field with respect to 
altitude, z, is of course assumed for each slab, z n , and associated with each such slab a set of four 
roots, £*, is found to exist. Two of the roots will exhibit a negative imaginary part (Im f nega- 
tive), corresponding to an upgoing propagation component {-\-z direction, fig. 2). Also, two 
of these roots will exhibit a positive imaginary part (Im f positive), corresponding to a down- 
going propagation component (—^-direction, fig. 3). Except for the topmost slab, it is neces- 
sary to consider both upgoing and downgoing components in this analysis. 

-> -> 

It is necessary [Johler, 1962] to equate the tangential E and H fields at each boundary, 

figure 2 (top and bottom of each slab), whereupon it can be concluded that a pXp + 2 matrix 
equation defines the propagation: 
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a l, 1 a l,2 a l, 3 a l, 4 a l. 5 a l, 6 
b l, 1 b l,2 b l, 3 b l,4 b l, 5 b l,6 
c l, 1 C l, 2 C l, 3 C l,4 C l, 5 c l, 6 
d l, 1 d l,2 d l,3 d l,4 d l,5 d l,6 

a 2, 3 a 2, 4 a 2, 5 a 2, 6 a 2, 7 a 2, 8 a 2, 9 a 2, 10 
b 2, 3 b 2,4 b 2,5 b 2,6 b 2,7 b 2. 8 b 2, 9 b 2, 10 
c 2, 3 c 2,4 c 2, 5 c 2, 6 C 2, 7 c 2, 8 °2, 9 C 2, 10 
d 2, 3 d 2,4 d 2 > 5 d 2 ) 6 d 2,7 d 2, 8 d 2, 9 d 2, 10 

a 3, 7 a 3, 8 a 3, 9 a 3, 10 a 3, 11 a 3, 12 a 3, 13 a 3, 14 
b 3, 7 b 3, 8 b 3, 9 b 3, 10 b 3, 11 b 3, 12 b 3, 13 b 3, 14 
c 3, 7 c 3, 8 C 3, 9 C 3, 10 c 3, 11 C 3, 12 C 3, 13 c 3, 14 
d 3,7 d 3, 8 d 3, 9 d 3, 10 d 3, 11 d 3, 12 d 3, 13 d 3, 14 

a p, p + 4 a p, p + 9 

P- P + 4 p, p + 9 

C P. P + 4 c p, p f 9 

P. P + 4 p, p + 9 
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(39) 



where the elements a n . . . can be deduced as a consequence of these boundary conditions? 
provided the reflection coefficients, T, and the transmission coefficients, U, T f are defined, 
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E/i 



(40) 



where, figure 3, n=l, 2, 3 . . ., p—1, p, and the four T f are defined as a transmission coeffi- 
cient into a topmost slab of infinite extent, z. 

Obviously, the matrix (39) can become quite large, and the success in handling such an 
equation (39) is dependent upon the memory capacity of a large-scale computer (such as the 
IBM-7090). For example, a 104X106 complex matrix equation (39) corresponds to the 
upper limit of 26 slabs used on the IBM-7090 and the CDC-1604 computers [Johler and 
Harper, 1962a and 1962b]. The size of such slabs were approximately 0.5 km. This was 
much finer detail than can ordinarily be resolved from N(h) curves of the lower ionosphere. 
Ordinarily, only 8 to 10 slabs exhibited quite remarkable convergence of the reflection coeffi- 
cients at frequencies of 10 to 100 kc/s for graphical significance in the computation. 

The distinction between the four propagation components, upgoing (i) and downgoing (r) , 
ordinary (o) and extraordinary (e), figure 3, becomes obscure when the ionosphere is con- 
sidered to be nonuniform in the vertical ^-direction. This is a consequence of the coupling at 

each boundary. Thus, the boundaries of each slab are introduced as an expression for the 

-> -> -> 

continuity of the tangential E and H fields and possibly the normal H at each boundary, 

figure 3, of the model plasmas. This is accomplished by equating the field immediately above 

and immediately below each boundary which, after considerable ado, results, for example, in 
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a 104X106 complex matrix equation, or 26 slabs [Jollier and Harper, 1962a and 1962b]. In 
addition to the transmission coefficients, U, four reflection coefficients referenced to the bottom 
of the lowest slab used (iV=10, h=65 in Houston composite model, figure 1, N=S, A = 45 in 
Blackout model), T ee , T em , T me , T mm can be defined. These coefficients describe the reflection 
into the region below the model plasma. Thus, T ee , used in this analysis, refers to vertical 
electric polarization of the incident wave and a corresponding vertical electric polarization of 
the reflected wave. T cm refers to the generation of the abnormal component (horizontal 
electric polarization) by the incident vertical electric wave. Similarly, T mm refers to vertical 
magnetic incident wave and vertical magnetic reflected wave, and T me refers to the correspond- 
ing abnormal component (horizontal magnetic). Employing this analysis procedure, and using 
the vertical electric reflection coefficient, T ee , three models were constructed as illustrated in 
figure 4. Model No. 1 employed the Houston Composite profile figure 2, and a distance of 
1609 km between transmitter and receiver with typical values of the other parameters. Model 
No. 2 is a local ionosonde blackout model appropriate to high latitude condition in the auroral 
zone. Model No. 3 is the same model as No. 1, but employs a greater distance beyond (he 
geometric-optical horizon. These models then define the propagation medium for the pulses 
to be described. 

3. Pulse Synthesis 

The phenomena of time separation of pulses into ground, first ionospheric reflection, etc., 
has been shown experimentally by Johler [1962] employing the radio-navigation-system pulse 
(Loran-C radio-navigation/timing system). Each such pulse can be synthesized as a super- 
position of three damped sinusoids as separate source-current moments. 



^(*)=£exp(-rf)- 



4 ex P (-^)— jexp (- 



•V2t), 



(41) 



where v=Ci-\-io> C} Pi=Ci+i(ci> c +2«j,), v 2 =Ci J ri(oo c — 2o) p ) and / c =w c /27r frequency, characteristic. 
By choice of various values for c x and oj p , a variety of radio-navigation type pulses can be 
investigated. Methods of inserting such pulses (damped sinusoids) in the Fourier transform- 
integrals and evaluating such integrals numerically and analytically have been described by 
Johler [1962] and Johler and Walters [1959]. Thus, both the damped sine wave and the damped 
cosine wave source current fields can be found b}^ evaluating the integral : 



E(t',d)=~f~\E(u,d)\ 

3 Olt'—<t>' c 



\ 



+i 



-f tan 



_! — (co c + co) 



Ci 



Vc?+(co c + c) 2 



cos — 

J + — — 



w£+<£c+tan 



-i — (o? c — co) T 



Vc?+(« c — co) 2 



sin 



[„«'-*:+tan-» ~ (( °; + a,) ] sin[- w <'+*;+tan-' - (M ;~ M) ] 



Vc?+(co c +co) 2 



\cl+(co c — oj) 2 



AJ 



c/oj (42) 



where 4>' c = — <t> c — 7r/2 = arg [E(u,d)\ and E(o),d)=E (co,d) or Ei(o* y d) for the groundwave or for 
the ionospheric wave respectively, where the field of the damped cosine source, Re F 8 (t)*is He 
E(t'jd) and the field of the damped sine source, — Im F s (t) is — Im E(t',d) or Re [iE(t'}d)\. 
The source-current moment exp ( — vt) is complex, since v=Ci-\-iti. This is analogous to 
the use of complex time function exp (iuf) in continuous time-harmonic wave analysis. The 
real part is implied, 



489 



s 

s 

3 I0" s 





1 1 1 1 






,. /^'^ 








yv \ 










A* \ 








/ / \ 












'/ 1 1 1 1 





I0 2 



I0 3 I0 4 

FREQUENCY, f, CYCLES /SECOND 



I0 5 



I0 6 



d ■ 1609 km 
<T " 5 
€ 2 = 80 

a ■ 0.85 

h * 65 km 

I ■ 60» 
H m = 0.5 
tl ■ 81.8° 
4> Q ■ 270* 
VERTICAL POLARIZATION 



MODEL NO 2 



d ■ 1416 km 

a = 5 

€ 2 = 80 

a = 0.85 

h » 45km 

I = 60 - 

H m - 0.5 

4>\ = 83.2* 

4>q ■ 270° 
VERTICAL POLARIZATION 



MODEL NO. 3 



d = 2125 km 

a ■ 0.005 

€ 2 ■ 15 

a ■ 0.85 

h = 65 km 

I ■ 67° 

H m ■ 0.505 

*o ■ 66« 

VERTICAL POLARIZATION 



— 



\ -25 



-50 — 



1 


1 I 1 




- 


\\\ 

MODEL NO. 2 V.A 

v- \ 


- 




V;\ 

t \ 

\\ \ 

MODEL NO 3 -\\ \ 












- 


k \ 

F- 1 
*• 1 

!'• 1 
1 • l 


- 


- 


i: 

i '• 
i : 
i : 


- 


1 


1 : 
1 : 
\ : 
1 : 
\ : 
l • 

1 1 1 


- 



I0 5 



I0 6 



490 



10 I0 2 I0 3 I0 4 

FREQUENCY, f, CYCLES /SECOND 

Figuee 4. Amplitude and phase of the complex transfer function, Ej(co,d), of the propagation medium for the 

first, j = l, ionospheric reflection. 

Three propagation medium models are illustrated: (1) The daytime-noon model for the lower geographic latitudes, distance = 1609 km; (2) 
The local ionosonde blackout model for higher geographic latitudes (auroral zone) distance = 1416 km; (3) The daytime-noon model for lower geo- 
graphic latitudes, distance =2125 km. 



Re exp (—vt) = exp (—Cit) cos (Sit) 

which is a damped sinusoid. The corresponding field, Re E(t' } d), is therefore implied in integral 
(42) as the oscillogram of the propagated pulse. The imaginary part of the integral, hence, 
corresponds to a source-dipole current moment, exp (— at) sin (Qt). Thus, in effect , two pulses 
are propagated simultaneously by the complex form of the source. This defines a source 
envelope, \F s (t)\, which for the case of sinusoids has a clear physical interpretation, (both for 
the field, \E(t',d)\, and the source, \F s (t)\) as the detected pulse envelope, corresponding to the 
results of an ideal detection process in a pulse receiver. Thus, the propagation of a third-type 
pulse wave or envelope wave is described by the integral (42). The form or shape of the en- 
velope wave is a precise description of the dispersion of the pulse, and methods of tagging a 
point-in-time on the pulse to be described below employ the quantity, \E(t' , d)\. It is of in teres! 
to note that the quantity, \E(t',d)\, and Re E(t',d) would be distinct even if the source, F s (t), 
were real. In other words, an envelope would be synthesized for the pulse by the integral 
(42) automatically. 

4. Tagging a Point-in-Time on the Pulse 

The measurement of the propagation time of a pulse or signal requires a method for tagging 
a point-in-time on the leading edge of the pulse. Sommerfeld [1914] and Brillouin [1914] 
established the consistency of signal velocity with relativity by defining such a point as the 
"time at which the signal takes on appreciable amplitude." Such a point-in-time while suffi- 
cient to prove the consistency of electromagnetic theory with the relativity principle, i.e., to 
prove that the signal always arrived at such times that the velocity of propagation is less than 
the speed of light, is inadequate for modern radio navigation and timing systems. Conse- 
quently, numerous precise experimental methods for defining and, indeed, tagging a point-in- 
time on a radio pulse have been devised by Johler [1963a], Doherty et al. [1901], and Frantz, 
Dean, and Frank [1957]. Although one method for tagging a point-in-time on the leading 
edge of a pulse can be more efficient experimentally than another, the method used mathe- 
matically is inconsequential as long as the mathematics does not lose computation precision 
in the tagging process. A point later than £'=0 is usually tagged since the signal amplitude is 
usually nil at zero time locally. The point tagged is inconsequential once it is tagged. Fur- 
thermore, the behavior of such a point-in-time on a pulse as distance, conductivity of the 
ionosphere or ground change can be studied. Such an analysis can lead to a criterion of pulse 
dispersion. 

The amplitude envelope minus the time derivative of the amplitude envelope method 
[Johler 1903a] was employed in this analysis: 

F{t\d)=Cmt'M~C 2 ^\E{t',d)\=--Q- (43) 

Thus, for given values of the constant, C 2 /Ci the root of the differential equation (43) was 
found, t' = T c , or the zero crossings of the function, F(t',d). The dispersion of the pulse can 
be studied in detail by examining several such points, C 2 /(\, as a function of the distance and 
conductivity parameters. 

5. Propagation Analysis of the Pulse 

The earliest pulse is illustrated, figure 5, at 1009 km (1000 statute miles); 1410 km (880 
statute miles); and 2125 km (1329 statute miles). These distances, especially the latter, were 
selected for checking the theory of propagation experimentally. Although the transmitters and 
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receivers for propagation over the distance, 2125 km, were located on land the intervening- 
surface was primarily sea water; hence, a conductivity, a=5, and a dielectric constant, e 2 =80, 
were employed for the groundwave, and a conductivity, cr=0.005, and dielectric constant, e 2 =15, 
were selected for the factor, F x . 

The groundwave pulses, Re E(f',d) reach a crest in approximately 5 cycles or 52 /xsec, 
at which point the field amplitude envelope, \E(t',d)\, is — 1.9(10~ 10 ) v/m, <i=2125 km, for the 
assumed dipole current moment, I l = l. This corresponds to a radiated power, P r , [Johler r 
1961] 

P,-1.6(10- 13 )co 2 (/ Z)/^w 

which at/ c =100 kc/s, z ~120tt, 

Pr — l.lilO'^) w peak radiated power. 

Obviously the field will increase linearly with dipole moments greater than unity or increase 
as to the square root of the power radiated, P T . The groundwave pulse used with the blackout 
model ionospheric wave was, of course, not affected by the condition of the ionosphere by 
definition of the groundwave pulse. The differences between the pulses are due only to the 
effect of distance and, hence, a greater or less effect of conductivity and earth curvature on the 
pulse. The vertical lapse factor, a, and the dielectric constant would only exhibit secondary 
effects; hence known approximate values for these parameters were assumed. 

The next pulse to arrive at the receiver, figure 6, after considerable ionospheric wave delay 
(20 to 40 /xsec) is the ionospheric wave pulse. The local time scale is now transformed to £'=£/, 
i.e., the earliest time the ionospheric pulse can arrive at the receiver over the ray, j= 1, figure 1, 
is ^'=0. The pulse corresponding to the 2125 km distance reaches the crest of the envelope 
in approximately 5 cycles or 58 /xsec as compared with 52 jusec for the corresponding ground- 
wave pulse, figures 5 and 6. 

The composite pulse can be formed since the pulse observed at a distant receiver is never 
a pure groundwave pulse or a pure ionospheric-wave pulse. There are areas of overlap [Johler, 
1962]. Writing the groundwave plus the ionospheric-wave pulse, 

UeE(t' } d)=Ke [JE (M)+#i (*',<*)], 

or using the same time scale, t', the pulse appears as illustrated in figure 7. The amplitude 
envelope \E(t',d)\ is also shown. The ionospheric-wave delay is shown. Note the short iono- 
spheric-wave 'delay (23.73 /xsec) for the blackout model. This is a consequence of the lowering 
of the ionosphere, figure 2. 

The amplitude envelope minus the time derivative of the amplitude envelope function, 
F(t* \d)y (34) is illustrated for the three pulses for the groundwave pulse, figure 8, for the iono- 
spheric-wave pulse, figure 9, and for the composite pulse, figure 10. The zero crossings were 
determined numerically for these pulses and the various points, C2/C1, tagged with precision 
are tabulated, tables 1, 2, 3, for the three models. The physical reality of the groundwave is 
obvious from a comparison of the pure groundwave tagged points-in-time 10.06 to 39.07 /xsec, 
table 1 for example, with the first tagged point on the composite pulse. Since these columns 
are identical there is no ionospheric wave contamination. However, when the pure ionospheric- 
wave pulse column is compared with the second tagged point-in-time on the composite wave, 
the contamination of the ionospheric-wave pulse by the ground- wave pulse can be noticed. 
The amount of signal delay for each tagged point-in-time on the ground and ionospheric-wave 
pulses is shown as a variable in the difference column as result of a change in the pulse dispersion. 

Table 2 corresponds to a high-geographic-latitude model which is characterized by the 
local ionosonde blackout during which ionization extends to lower levels. The first tagged 
point-in-time becomes contaminated quite early on the leading edge of the composite pulse 
with the ionospheric-wave pulse. Thus, the point, 24.57 /xsec includes some of the ionospheric 
wave pulse. It is necessary to back up to 16.19 /xsec to obtain a pure ground wave. It should 
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Table 1. Model 1: low-latitude, daytime-noon ionosphere 
d=1609km (1000 statute miles) 





<r=5mhos/m 
e 2 =80 


/ e =100kc/s 
7=60° 


Ci=3 (10)4 
/ P =2.5kc/s 


0a = 27O° 

/i=65km 
















did 




Groundwave 


Ionospheric 
wave 


Ionospheric 
wave 41.16 a 


Composite ground-ionospheric wave 
tagged point-in-time 




First 


Second 


Difference 


5(10-6). 


10.06 
16.35 
25.11 
30.69 
34.34 
37.04 
39.07 


t,nsec 
52.13 
59.02 
67.88 
73.41 
77.15 
79.84 
81.86 


usee 
10.97 
17.86 
26.72 
32.25 
35. 99 
38. G8 
40.70 


10.06 
16.35 
25.11 
30.69 
34.34 
37.04 
39.07 






10-5 


60.67 
69.48 
74.87 
78.54 
81.16 
83.15 


44.32 


2(10-5) 


44.37 


3(10-5) 


44.18 


4(10-5) 


44.20 


5(10-5) 


44.12 


6(10-5) 


44.08 







1 Ionospheric wave delay relative to the groundwave— 41.16 Msec; continuous wave: 

Amplitude Phase 

Ground 6.6213(10-9) v/m -2.1255 radians 

Ionospheric 6.1096(10-8) v/m 1.7517 radians 



able 2. Model 2: high-latitude disturbed ionosphere 
d = 1416 km (880 statute miles) 



o-=5mhos/m 

€2 = 80 



/=100kc/s 
/=60° 



ci =3(104) 

/ P =2.5kc/s 



<t> a =270° 
h = 45 km 



CilCx 



Groundwave 



Ionospheric 
wave 



Ionospheric 
wave-23.73 a 



Composite ground-ionospheric wave 
tagged point-in-time 



First 



Second 



Difference 



10-5.... 
2(10-5). 
3(10-5) . 
4(10-5). 
5(10-5). 
6(10-5). 



16. 19 
24.98 
30.52 
34.21 
36.89 
38.92 



t, usee 
46.00 
54.84 
60.36 
64.09 
66.77 
68.79 



usee 
22.27 
31.11 
36.63 
40.36 
43.04 
45.06 



16.19 
24.57 
29.28 
30.61 
31.04 
31.40 



58.97 
64.45 
68.09 
70.71 
72.67 



34.40 
35.17 
37.48 
39.67 
41.27 



a Ionospheric wave delay relative to the groundwave— 23.73 /usee; continuous wave: 
Amplitude Phase 

Ground 1.1060(10-8) v/m -1.8410 radians 

Ionosphere 3.7966(10-8) v/m -2.8235 radians 



be noted that the composite wave after the arrival of the ionospheric wave contains some ground 
wave. This contamination with ground wave has been observed experimentally by Doherty 
[1963] in the Aleutian Islands and in the North Atlantic at distances between 600 and 900 
statute miles employing the Loran-C radio navigation/timing system. Since the ionospheric- 
wave delay is, in general, less in the arctic regions or northern latitudes, a reflection process 
at a lower altitude is indicated. The consequence of a low-level ionosphere reflection as far 
as the automatic equipment is concerned is illustrated by figure 10. Note the tagged point 
for the Blackout Model (10~ 5 ), for example. The null-seeking devices which search for the 
zero crossing of the amplitude envelope minus the time derivative of the amplitude envelope 
function would never find such a point on the ionospheric pulse unless the operator reset 
C 2 /Ci to 2(10~ 5 ). Indeed, the equipment would, if an automatic search were permitted, lock 
onto the ground-wave near t' =16 jusec during such a blackout. 

Table 3 illustrates the low-geographic latitude quiescent ionosphere at great distance of 
2125 km which is beyond the geometric-optical horizon, figure 1. Here the amount of ground- 
wave contamination of the second point tagged on the composite pulse becomes smaller as a 
result of the high attenuation of the groundwave pulse, if a point is selected quite high on 
the pulse, <7 2 /Ci = 6(10- 6 ) or £' = 80.82 M sec. 
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Figure 11. The zero crossings of the cycles on the ionospheric-wave, ground-wave, and composite pulse, illustrating 
cycle distortion or apparent cycle time shift on the composite pulse, Models 1, 2, 3. 
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Table 3. Model 3; low-latitude, daytime-noon ionosphere 

d=2125 km (1321 statute miles) 

/ p =2.5 kc/s 7=67° 
€1 = 80 4> a = 66° 

C2 = 15 ft = 65km 



f e = 100 kc/s 
Ground wave 
Ionospheric wave 



ci =3(104) 
<7=5mhos/m, 
o- = 0.005 mhos/m, 



C2/C1 


Groundwave 


Ionospheric 
wave 


Ionospheric 
wave— 41.13 a 


Composite ground-ionospheric wave 
tagged point-in-time 




First 


Second 


Difference 


5(10-6) 


10.58 
16.73 
25.55 
31.12 
34.79 
37.41 
39.45 


t, usee 
56.00 
62.22 
71.17 
76.67 
80.46 
83.11 
85.12 


usee 
14.87 
21.09 
30. 04 
35.54 
39.33 
41.98 
43.99 


10.58 
16.73 
25.55 
31.12 
34.79 
37.40 
39.70 






10-5 






2(10-5) 


64.60 
71.39 
75.64 
78.62 
80.82 


39.05 


3(10-5) 


40.27 


4(10-5) 


40.85 


5(10-5) 


41.22 


6(10-5) 


41.12 







1 Ionospheric wave delay relative to the groundwave— 41.13; Msec continuous wave: 
Amplitude Phase 

Ground 1.7311(10-9) v/m -2.8871 radians 



Ionospheric 7.5309(10-8) v/m 



-1.1389 radians 



The corresponding time-harmonic wave solution for the ground and ionospheric wave is 
shown at the bottom of each table (1, 2, and 3). The phase lag of the groundwave on table 
1 in microseconds is [|-2.1255|/6.28318]10 = 3.382 /*sec for the groundwave to which 10 /*sec 
or 2w radians can be added, getting 13.382; 23.382 . . . /*sec. This phase is ambiguous by 
±27r without the pulse detection procedure outlined in this paper, i.e., without a precise method 
for tagging a pulse in the time domain. The phase of the continuous time-harmonic wave 
becomes even more obscure when applied to the composite ground and ionospheric wave. 
Finally the zero crossings of the individual cycles under the envelope of the pulse are illustrated, 
figure 11. The curves are approximately straight lines with a 45° slope. Indeed, this would 
be the case without cycle distortion. However, the composite wave exhibits some cycle 
distortion in the vicinity of overlap of the two pulses. Also, at early times on the pulse, the 
first cycle is distorted appreciably as shown in the 2125 km example, for both the ground and 
the ionospheric wave pulses. 

6. Conclusions 

The methods which employ the direct analytic or numerical procedures for the evaluation 
of the Fourier transform-integral for the groundwave pulse can be applied to ionospheric 
wave pulses. The form or shape of the composite pulse compares very well with those ob- 
served in practice. The phenomenon of time separation of pulses observed experimentally 
has been demonstrated mathematically. The complete detailed analysis of pulse dispersion, 
signal velocity, and ionospheric wave pulse or groundwave pulse timing requires precise methods 
for tagging a point-in-time on the pulse — both mathematical and experimental. One such 
mathematical method has been successfully applied in this analysis. 
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